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ABSTRACT 

If the origin of life and the evolution of observers on a planet is favoured by atypical properties of 
a planet's host star, we would expect our Sun to be atypical with respect to such properties. The 
Sun has been described by previous studies as both typical and atypical. In an effort to reduce this 
ambiguity and quantify how typical the Sun is, we identify 11 maximally-independent properties that 
have plausible correlations with habitability, and that have been observed by, or can be derived from, 
sufficiently large, currently available and representative stellar surveys. By comparing solar values for 
the 11 properties, to the resultant stellar distributions, we make the most comprehensive comparison 
of the Sun to other stars. The two most atypical properties of the Sun are its mass and orbit. The Sun 
is more massive than 95 ± 2% of nearby stars and its orbit around the Galaxy is less eccentric than 
93 ± 1% of FGK stars within 40 parsecs. Despite these apparently atypical properties, a x^-analysis of 
the Sun's values for 11 properties, taken together, yields a solar Xq = 8.39 ± 0.96. If a star is chosen 
at random, the probability that it will have a lower value (~ be more typical) than the Sun, with 
respect to the 11 properties analyzed here, is only 29± 11%. These values quantify, and are consistent 
with, the idea that the Sun is a typical star. If we have sampled all reasonable properties associated 
with habitability, our result suggests that there are no special requirements for a star to host a planet 
with life. 

Subject headings: stars: fundamental parameters — stars: statistics — Sun: fundamental parameters 
— Sun: general 
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1. INTRODUCTION 

If the properties of the Sun are consistent with the idea 
that the Sun was randomly selected from all stars, this 
would indicate that life needs nothing special from its 
host star and would support the idea that life may be 
common in the universe. More particularly, if there is 
nothing special about the Sun, we have little reason to 
limit our life-hunting efforts to planets orbiting Sun-like 
stars. As an example of the type of anthropic reasoning 
we are using, consider the following situation. Suppose 
uranium (a low abundance element in the Solar System 
and in the universe) was central to the biochemistry of 
life on Earth. Further, suppose that a comparison of 
our Sun to other stars showed that the Sun had more 
uranium than any other star. How should we interpret 
this fact? The most reasonable way to proceed would 
be to try to evaluate the probability that such a coin- 
cidence happened by chance and to determine whether 
we are justified in reading some importance into it. Al- 
though a correlation does not necessarily imply cause, 
we think that a correlation between the Sun's anoma- 
lous feature and life's fundamental chemistry would be 
giving us important clues about the conditions necessary 
for life. Specifically, the search for life around other stars 
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as envisioned by the NASA's Terrestrial Planet Finder or 
ESA's Darwin Project and as currently underway with 
SETI would change the strategy to focus on the most 
uranium-rich stars. Another example: Suppose the Sun 
had the highest [Fe/H] of all the stars that had ever been 
observed. Then high [Fc/H] would be strongly implicated 
as a precondition for our existence, possibly by playing a 
crucial role in terrestrial planet formation. These are ex- 
aggerated examples of the more subtle correlations that a 
detailed and comprehensive comparison of the Sun with 
other stars could reveal. 

Whether the Sun is a typical or atypical star with re- 
spect to 1 or a few properties has been addressed in previ- 
ous studies. Using an approach similar to ours (compar- 
ing solar to stellar properties from particular samples), 
some studies have suggeste d that th e Sun is a typical 
star ( jCustafsson 1998 ; AUende Prietd ,2006). while other 
studies ha ve sug gested that the Sun is an atypical star 
(|Gonzaledll999allbl : lGonzalez et al.ll2001h . This apparent 
disagreement arises from three problems: 

1. Language used to describe whether the Sun is, or 

is not typical, is often confusingly qualitative. For 
example, reporting the Sun as "metal-rich" , can 
mean that the Sun is significantly more metal-rich 
than other stars (e.g. more metal-rich than 80% of 
other stars) or it can mean that the Sun is insignif- 
icantly metal-rich (e.g. more metal-rich than 51% 
of other stars). 

2. Selection effects: Stellar samples chosen for the com- 

parison can be biased with respect to the property 
of interest. 



3. Inclusion (or exclusion) of stellar properties for which 
it is suspected or known that the Sun is atypical, 
will make the Sun appear more atypical (or typi- 
cal). 

In this paper we address problem 1 by using only quan- 
titative measures when comparing the Sun's properties 
to other stars. Our main interest is to move beyond 
the qualitative assessment of the Sun as either typical 
or atypical, and obtain a more precise quantification of 
the degree of the Sun's (a)typicality. In other words, we 
want to answer the question 'How typical is the Sun?' 
rather than 'Is the Sun typical or not?' There are at 
least two ways to quantify how typical the Sun is. This 
can be done for individual parameters by determining 
how many stars have values below or above the solar 
value (Table [3]). This can also be done by a joint analy- 
sis of multiple parameters (Table [2]) ■ If there are several 
subtle factors that have some influence over habitabil- 
ity, a quantitative joint analysis of the Sun's properties 
may allow us to identify these factors without invoking 
largely speculative arguments linking specific properties 
to habitability. 

With respect to problem 2, most previous analyses 
have compared the Sun to subsets of Sun-like stars se- 
lected to be Sun-like with respect to 1 or more param- 
eters. In such analyses, the Sun will appear typical 
with respect to any parameter(s) correlated with 1 of 
the pre-selected Sun-like parameters. For example, el- 
emental abundances [X/H] are correlat ed with metal- 
hcityS [Fe/H]. The sample of Edvardsso n et all (|1993aD 
was selected to have a wide range of [Fe/H]. This pro- 
duced a metallicity distributio n unrepresentative of stars 
in general. Recognizing this, lEdvardsson et al.1 (|1993aD 
conditioned on solar metallicity, [Fe/H] w and then 
compared solar abundances for 12 elements to the abun- 
dances in a group of nearby stars with solar iron abun- 
dance, solar age and solar galactocentric radius. They 
found the Sun to be "a quite typical s tar for its inetallic - 
ity, age and galactic orbit" . Similarlv. lGustafssonI (|1998f ). 
after comparing various properties of the Sun to solar- 
type stars (stars of similar mass and age) , concluded that 
the Sun seems very normal for its mass and age; "The 
Sun, to a remarkable degree, is solar type" . The stellar 
samples we use for comparison with the Sun are, in our 
judgement, the least-biased samples currently available 
for such a comparison. 

To address problem 3, in Section[2]we compare the Sun 
to other stars using a large number (11) of maximally- 
independent properties with plausible correlations with 
habitability. These properties can be observed or derived 
for a sufficiently large, representative stellar sample (Ta- 
ble[T]). Any property of the Sun or its environment which 
must be special to allow habitability would show up in 
our analysis. However, in contrast to previous analyses 
which have looked for solar anomalies with respect to 
individual properties, we perform a joint analysis that 
enables us to quantify how typical the solar values are, 
taken as a group. In Section [3l the differences between 
the solar values and the stellar samples' medians are used 
to perform first a simple and then an improved version 

^ Metallicity: [Fe/H] is the fractional abundance of Fe rel- 



of a x^-analysis to estimate whether the solar values are 
characteristic of a star selected at random from the stellar 
samples. The results of our joint analysis are presented 
in Figure [T51 of Section 2] We find that the solar values, 
taken as a group, are consistent with the Sun being a ran- 
dom star. However, there are important caveats to this 
interpretation associated with the compromise between 
the number of properties analyzed, and their plausibil- 
ity of being correlated with habitability. In Sections \5\ 
and [6] we discuss these caveats and summarize. We dis- 
cuss the levels of correlation between our 11 properties 
in Appendix A. 

2. STELLAR SAMPLES AND SOLAR VALUES 

We are looking for a signal associated with a prerequi- 
site for, or a p roperty that favor s, the origin and evolu- 
tion of hfe fsee lGustafssonllT998l for a brief discussion of 
this idea) . If we indiscriminately include many properties 
with little or no plausible correlation with habitability, 
we run the risk of diluting any potential signal. If we 
choose only a few properties based on previous knowl- 
edge that the Sun is anomalous with respect to those 
properties, we are making a useful quantification but we 
are unable to address problem Hi. We choose a middle 
ground and try to identify as many properties as we can 
that have some plausible association with habitability. 
This strategy is most sensitive if a few unknown stellar 
properties (among the ones being tested) contribute to 
the habitability of a terrestrial planet in orbit around a 
star. 

An optimal quantitative comparison of the Sun to 
other stars would require an unbiased, large represen- 
tative stellar sample from which independent distribu- 
tions, for as many properties as desired, could be com- 
pared. Such a distribution for each property of interest 
would allow a straightforward analysis and outcome: the 
Sun is within the n% of stars around the centroid of the 
TV-dimensional distribution. However, observational and 
sample selection effects prevent the assembly of such an 
ideal stellar sample. 

In this study, we compare the Sun to other stars with 
respect to the following 11 basic physical properties: 
(1) mass, (2) age, (3) metaUicity [Fe/H], (4) carbon- 
to-oxygen ratio [C/O], (5) magnesium-to-sihcon ratio 
[Mg/Si], (6) rotational velocity wsini, (7) eccentricity 
of the star's galactic orbit e, (8) maximum height to 
which the star rises above the galactic plane Z^nax, (9) 
mean galactocentric radius Rcai, (10) the mass of the 
star's host galaxy Mgai, (11) the mass of the star's host 
group of galaxies Mgi-oup- These 11 properties span a 
wide range of stellar and galactic factors that may be as- 
sociated with habitability. We briefly discuss how each 
parameter might have a plausible correlation with hab- 
itability. For each property we have tried to assemble 
a large, representative sample of stars whose selection 
criteria is minimally biased with respect to that prop- 
erty. For each property the percentage of stars with 
values lower and higher than the solar value are com- 
puted. For properties (9), (10) and (11), the uncertain- 
ties in the percentages are determined from the uncer- 
tainties of the distributions. For the rest of the proper- 
ties, nominal uncertainties A, on the percentag es were 
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TABLE 1 

Samples used to produce the stellar distributions plotted in Figi 
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Uliaracteristic width of distribution in the direction of the solar value. 



Property 


Range 


Medi< 

Ml/2 


Mass [Mq] 


0.08 - 2 


0.:. 


Age [Gyr] 


0-15 


5. 


[Fe/H] 


-1.20 - +0.46 


-0. 


[C/O] 


-0.22 - +0.32 


o.c 


[Mg/Si] 


-0.18 - +0.14 


o.c 


V sin i [km s -*-] 


0-36 


2.t 


e 


0-1 


0.1 


.Zmax [kpc] 


- 9.60 


0.1 


R-Gal [kpc] 


0-30 


4. 


Mgal [Mq]^ 


10^ - 10^2 


lOi 


Mgroup [Mq]'' 


IqS - 1013 


lOi 



t 

''IWright et al.' (WO^), ( see footnote in Sec. I2.2| |. 

^ G99: Gustafsson et al. 19ji, R03: [Re ddv et at pOOl . BF06: IBensbv fc Feltzinell200g . 
R03: Redd y et al. .2003. BOS: IBensW et al.n2005[ 

Solar rotational velocity corrected for random inclination (see Sec. I2.5| l. 

Sub-set of stars within the mass range: 0.9 Mq < M < 1.1 Mq. 
' Calculated using the solar galactic motion lIDehnen fc Binnevlll99^ and the Galactic potential (see Sec. I2.6| |. 
s Sub-set of volume complete A5-K2 stars within 40 pc. 

^ Integrated solar orbit in the Galactic potential of lFlynn et al.l II1996I ) (see Sec. I2.6l l. 

' [Eisenhaucr et al. 200j; 

^ BS80: Bahcall fc S oneir a| [T980l . G96: IGould et aL|[T99l E05: lEisenhauer et al.|[2005l . 
^ Stellar mass, not total baryonic mass, nor total mass. 

' D94: [Driver et ahlirgp . CB99: ICourteau fc van den Bergb|[T999l, LOO: ILovedavll2000l . BJOl: iBelT 



: de Jon a i200ll . J03: IJarrett et an[200a . 



the fraction of stars with a lower (higher) value than the 
Sun and A'tot is the total number of stars in the sample. 
The solar value is indicated with the symbol "0" in all 
figures. 

We compare the Sun and its environment to other stars 
and their environments. The analysis of these larger en- 
vironmental contexts provides information about proper- 
ties that otherwise could not be directly measured. For 
example, suppose the metallicity of the Sun were normal 
with respect to stars in the solar neighborhood but that 
these stars as a group, had an anomalously high metal- 
licity with respect to the average metallicity of stars in 
the Universe. This fact would strongly suggest that hab- 
itability is associated with high metallicity, but our com- 
parison with only local stars would not pick this up. In 
the absence of an [Fe/H] distribution for all stars in the 
Universe, we use galactic mass as a convenient proxy for 
any such property that correlates with galaxy mass. 

2.1. Mass 

Mass is probably the single most important charac- 
teristic of a star. For a main sequence star, mass 
determines luminosity, effective temperature, main se- 
quence life-time and the dimensions, UV insolation and 
temporal stability o f the circumstellar habitable zone 
(|Kasting et al.lll993D . 

Low mass stars are intrinsically dim. Thus a com- 
plete sample of stars can only be obtained out to a dis- 
tance of ^ 7 parsecs (« 23 lightyears). Figure [1] com- 
pares the mass of the Sun to the stellar mass distribu- 
tion of the 125 nearest main sequence stars within 7.1 pc , 
as compiled by the RECONS consortium (Henry '2006) . 
Over-plotted is the stellar Ini tial Mass Fun ction (IMF; 
see eqs. 4 & 5 and Table 1 of lKroupall2002f l normalised 




Fig. 1. — Mass histogram of the 125 nearest stars l|Henrvll2006l . 
RECONS). The median (/^i/j = 0.33 Mq) of the distribution is 
indicated by the vertical grey line. The 68% and 95% bands around 
the median are indicated respectively by the vertical dark grey and 
light grey bands. We also use these conventions in Figs. 2-11. The 
solid curve and hashed area around it represents the Initial Mass 
Function (IMF) and its associated uncertainty (Kroupa 2002). The 
Sun, indicated by "0", is more massive than 95 ±2% of these stars. 



sal ()Kroupa fc Weidnerll2005f) . these nearby comparison 
stars are representative of a much larger sample of stars. 
There is good agreement between the histogram and the 
IMF — the Sun is more massive than 95 ± 2% of the 
neares t stars, and mo re massive than 94 ±2% of the stars 
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sample. Including them yields 94% — the same result 
obtained from the IMF. Our 95% ± 2 % result should b e 
compared with the 91% reported by iGonzale j (|1999bf ). 
The Sun's mass is the most anomalous of the properties 
studied here. 

2.2. Age 

If the evolution of observers like ourselves takes on av- 
erage many billions of years, w e might expect the Sun 
to be anomalously old |Carteii ri983). Accurate estima- 
tion of stellar ages is difficult. For large stellar surveys 
(> a few hundred stars), the most commonly used age 
indicators are based on isoch rone fitting and/or c hromo - 
spheric activity (i?HK index). iRocha-Pinto et all (|2000bf ) 
have estimated a Star Formation Rate (SFR), or equiv- 
alently, an age distribution for the local Galactic disk 
from chromospheric ages of 552 late-type (F8-K2) dwarf 
stars in the mass ran ge 0.8 M(t, < < 1.4 Mq at 



distances d < 200 pc (|Rocha-Pinto et al.ll2000ah . They 
applied scale-height corrections, stellar evolution correc- 
tions and volume incompleteness corrections that con- 
verted the observed age distribution in to the total num- 
ber of star s born at any g i ven t ime. [Hernandez et all 
(|2000D and iBertelli fc Nasil (|200in have made estimates 
of the star formation rate in the solar neighborhood 
and favou r a smoother dis tribution (fewer bursts) than 
[Rocha-Pi nto et all (|2000bD . 

In Figure [5] we compare the chromospheric age of the 
Sun (tq = 4.9 ± 3.0 Gyr, iWright et al.][200l ^ to the 
stellar age distribution re presenting the Galactic SFR 
(jRocha-Pinto et al.l[2000bh . The median of this distri- 
bution is 5.4 Gyr. The Sun is younger than 53 ± 2% 
of the stars in the thin disk of our Galaxy. Over-plotte d 
is the cosmic SFR derived bv lHopkins fc BeacomI |2006D . 
According to this distribution with a median Hi/2 = 9.15 
Gyr, the Sun was born after 86 ± 5% of the stars that 
have ever been born. 

The Galactic and cosmic SFRs are different because 
the cosmic SFR was dominated by bulges and elliptical 
galaxies in which the largest fraction of stellar mass in the 
Universe resides. Bulges and elliptical galaxies (early- 
type galaxies) formed their stars early and quickly and 
then ran out of gas. The disks of spiral galaxies, like 
our Milky Way, seem to have undergone irregular bursts 
of star formation over a longer period of time as they 
interacted with their satellite galaxies. 

The volume li mited ((i,„ax = 40 pc) sub-set from 
iNordstrom et al.l (j2004t ) contains isochrone ages for 1126 
A5-K2 stars. The median of this sub-set is 5.9 Gyr and 
the Sun is younger than 55 ± 2% of the stars. The sim- 
ilarity of this isochrone age result to the chromospheric 
age result is not obvious since the agreement between 
these two age technique s is rather poor. T his mismatch 
ca n be seen in Fig.fT^ . lReid et"all ()2007f ) . and in Fig. 8 
of lFeltzing et all (| 20011 ). 

2.3. Metallicity 

Iron is the most frequently measured element in nearby 
stars. Metallicity [Fe/H], is known to be a proxy for 

''' To ensure that the Sun's age is determined in the same way 
as the stellar ages to which it is being compared, we adopt the 
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Fig. 2. — The Galactic stellar age distribution (median /^i/2 = 
5.4 Gyr) from Ro cha-Pinto et al.l I l2000bl '). The Sun is younger than 
53 it 2% of the stars in the disk of our Galaxy. The grey curve is the 
cosmic Star F ormation Rate (SFR) with its associated uncertainty 
IIHopkins &: Bcacom 2006). according to which the Sun is younger 
than 86 it 5% of the stars in the Universe. 



the fraction of a star's mass that is not hydrogen or 
helium. In the Sun and possibly in the Universe, the 
dominant contributors to this mass fraction in order 
of abundance are: 0(44%), C(18 %), Fe(10%), Ne(8%) 
Si(6%), Mg(5%) N(5 %), S(3%) (lAsplund et all 120051 : 
iTruran fc Hegej|2005f ). The corresponding abundances 
by number are: 0(48%), C(26%), Ne(7%), N(6%), Mg 
(4%), Si(4%), Fe(3%), S(2%). Importantly for this anal- 
ysis, this short list contains the dominant elements in the 
composition of terrestrial planets (O, Fe, Si and Mg) and 
life (C, O, N and S). 

Over the last few decades, much effort has gone into 
determining abundances in nearby stars for a wide range 
of elements. Stellar elemental abundances for element 
X are usually normalised to the solar abundance of 
the same element using a logarithmic abundance scale: 
[X/H]* = log(X/H)^ - log(X/H)0. Hence all solar ele- 
mental abundances [X/H]q, are defined as zero. Spec- 
troscopic abundance analyses are usually made differen- 
tial relative to the Sun by analysing the solar spectrum 
(reflected by the Moon, asteroids or the telescope dome) 
in the same way as the spectrum of other stars. In this 
approach, biases introduced by the assumption of local 
thermodynami c equilibrium (LTE), larg ely cancel out for 
Sun-hke stars (jEdvardsson et al.]|1993bl ). 

A comparison between solar and stellar iron abun- 
dances is a common feature of most abundance sur- 
veys and most have conclud ed that the Sun is metal- 
rich co mpared to other stars (|Gustafssonlll998l : iGonzalej 
Il999allbl ). However, for our purposes, the appropriateness 
of these comparisons depends on the selection criteria of 
the stellar sample to which the Sun h as been compared. 
Stellar metalli city an a lyses s uch as lEdvards son et ahl 
(1993 a): , Redd v et all (l2003f) : INordstrom eTal, (2004f) : 
iValenti fc Fisched ^200^ have stellar samples selected 
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[Fe/H] 

Fig. 3. — Stellar met allicity histogr am of the 453 FGK Hippar- 
cos stars within 25 pc (IGrether & Lin eweaver 2007i). The median 
IJ.l/2 = —0.08. The Sun is more metal-rich than 65 ± 2% of the 
stars. 



the Galaxy and their sample is biased towards low 
metallicity (average [Fe /Hl= —0.25). The sample of 
IValenti fc Fischeil l2005l (average [Fe/H]= -0.01), was 
selected as a planet candidate l ist and contains some 
bias t owards high metallicity (see IGrether fc Lineweaverl 
[200l To assess how typical the Sun is, Gustafssoni 
()1998f l limited the sample of Edvardsson et al. (1993a) 
to stars with galactocentric radii within 0.5 kpc of the 
solar galactocentric radius, and to ages between 4 and 6 
Gyrs. The distribution of stars given by this criteria has 
an averag e [Fe/H]= —0.0 9. 

Grcther fc Lineweaveil (|2006l l2007f ) compiled a sample 
of 453 Sun-like stars within 25 pc. These stars were se- 
lected from the Hipparcos catalogue, which is essentially 
complete to 25 pc for stars within the spectra l type rang e 
F7~K3 and absolute magnitude of My = 8.5 (iRei J2002l ) . 
Metallicities for this sample were assembled from a wide 
range of spectroscopic and photo metric surveys. In Fig- 
ure [21 we compare the Sun to the IGrether fc Lineweaveil 
(jiOO?) sample, which has a median [Fe/H]= -0.08. To 
our knowledge this is the most complete and least-biased 
stellar spectroscopic metallicity distribution. The Sun is 
more metal-rich than 65 ± 2% of these stars^ 

Thi s result should be compared with iFavata et al.l 
()1997( ) who constructed a volume-limited (dmax = 25 
pc) sample of 91 G and K d warfs ranging in col or index 
{B ~ V) between 0.5 - 0.8 (jFavata et al.lll996f ). Their 
distribution has a median [Fe/H]= —0.05 and compared 
to this sa mple, the Sun is m ore metal rich than 56 ±5% of 
the stars. iFuhrmaiml ()2008D compared the Sun to a vol- 
ume complete (dmax = 25 pc) sample of about 185 thin- 
disk mid-F-type to early K-type stars down to My = 6.0. 
He finds a mean [Fe/H] = -0.02±0.18. This mean [Fe/H] 
is lowered by 0.01 dex if the 43 double- lined spectroscopic 
binaries in his sample are included. His results are con- 
sistent with ours. 



The elemental abundance ratios of a host star have 
a major impact on its proto-planetary disk chemistry 
and the chemical compositions of its planets. Oxy- 
gen and carbon make up ~ 62% of the solar system's 
non- hydrogen- non- h elium mass content (Z ~ 0.0122, 
lAsplund et al.ll2005l ). Carbon and oxygen abundances 
are among the hardest to determine. This is due to 
high temperature sensitivity and non-LTE effects in their 
permitted hues (e.g. C I A6588, O I A7773), and to 
the presence of blends in the forbidden lii ies ([ C I] 
A8727, [O I] A6300) . SeelAllende Prieto et all (I2OO1I ) and 
iBensbv fc Feltzinei ^200^ for details on C and O abun- 
dance derivations. 

Carbon pairs up with oxygen to form carbon monox- 
ide. In stars with a C/0 ratio larger than 1, most of the 
oxygen condenses into CO which is largely driven out of 
the incipient circumstellar habitable zone by the stellar 
wind. In this oxygen-depleted scenario, planets formed 
within the snow line are formed in reducing environments 
and are mostly compo sed of carbon compounds , for ex- 
ample silicon carbide (jKuchner fc Seage"i1l2005| ). Thus, 
the C/0 ratio could be strongly associated with habit- 
ability. 

As most heavy element abundances relative to hydro- 
gen (e.g. [0/H], [C/H], [N/H]) are correlated with [Fe/H], 
they were not included in our analysis. After the over- 
all level of metallicity (represented by [Fe/H]), and after 
the ratio of the two most abundant metals, [C/0], the 
magnesium to silicon ratio [Mg/Si] is the most important 
ratio of the next most abundant elements (excluding the 
noble gas Ne). For example [Mg/Si] sets the ratio of 
olivine to pyroxene which determines the ability of a sil- 
icate mantle to retain water (H. O'Neill 2007, private 
communication) . 

Stellar elemental abundance ratios are defined as 
[X1/X2]* — [Xl/ll]^, — [X2/H]^. Hence, systematic er- 
rors associated with the determination of absolute solar 
abundances cancel for abundances relative-to-solar. We 
compile [C/0] and [Mg/Si] ratios from samples with the 
largest number of stars and highest signal-to-noise ratio 
stellar spectra: 

[C/ 0].— 256 stars from iGustafsson et all 119991 : 
iReddv et al.ll2003l : iBensbv fc Feltzing||2006l 

[Mg /Si].— 231 stars from IReddv et al.l 12003) : 
IBensbv et al.ll2005l 

Due to their selection criteria, these samples are biased 
towards low metallicity and therefore cannot be used to 
create a representative [Fc/H] distribution. Because a 
correlation exist s between the [C/0] an d [Mg/Si] ratios 
and [Fe/H] (e.g. IGustafsson et al.lll999l ). the samples we 
use have a relatively narrow range of [Fe/H] to reduce 
the influence of the correlation. Therefore, these small 
correlations can be neglected in this study — see the bot- 
tom panels of Fig. |4| where [Fe/H] versus [C/0] as well 
as [Fe/H] versus [Mg/Si] are plotted. The top panels 
show the corresponding stellar distribution histograms. 
The Sun's [C/0] ratio is lower than 81 ib 3% of the stars. 
This is consistent with Gonzalez (19 99br) w ho suggested 
— based on data from Edvardsson et al.l (|1993ai) and 
IGustafsson et all ()1999D — that the Sun has a low [C/0] 
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4. — A: Comparison of the Sun's carbon-to-oxygen ratio ([C/0]q = 0) to the [C/O] ratios of 256 stars compiled from lGustafsson et ahl 



Dmparison 

||T999); Red dv et all 1(2003 ) and IBensb-TfcF eltzing (2006). The Sun's [C/O] ratio is lower than 81 ± 3% of the stars in this sample which 
has a median /ii/2 = 0.07. B: Comparison of the Sun's magnesium-to-silicon ratio ([Mg/Si]© = 0), to [Mg/Si] values from 231 stars from 
I Reddv et al.l l|2003l ) and IBensbv et aP l|2005l ) . The Sun's [Mg/Si] ratio is lower than 66 ± 3% of the stars in this sample with median 
^1/2 = 0.01. The bottom panels C & D show the small correlations of these distributions with [Fe/H]. These small correlations can be 
neglected for this study. 



the Sun is oxygen poor compared to solar metallicity 
stars. 

The Sun's [IVIg/Si] ratio is lower than 66 ± 3% of the 
stars. The [C/O] and [IVIg/Si] ratios are also largely in- 
dependent of each other (see Fig.[T3]in Appendix A). 

2.5. Rotational Velocity 

Stellar rotational velocities are related to the specific 
angular momentum of a protoplanetary disk and possibly 
to the magnetic field strength of the star during planet 
formation, and to protoplanetary disk turbulence and 
mixing. An unusually low stellar rotational v elocity may 
be associated with the presence of planets (jSoderbloml 
[1983,). One or several of these factors could be related to 
habitability. 



p. 485). In order to minimise the effect of this correla- 
tion (and maximize independence between parameters), 
we assembled a sample containing 27 6 stars within the 
mass range 0.9-1.1 Mq (F8-K2) from 
(2005"). The selection criteria of the 



Valenti fc Fischer 



Valenti fc Fischer 



2005 ) stars introduces some bias against more ac- 
tive stars. We compa red the high vsmi tail of our 
IValenti fc Fische^ (l20(35|) sample with th e high v sin i tail 
of a sub-sample fr om iNordstrom et al.l (120 04') . We esti- 
mate that for our IValenti fc Fischer] (|20()5 ) sample, the 
bias introduced by the selecti on criteria is lower than 
~ 5%. The vsini values in IValenti fc Fischer! (|2005f ) 
are obtained by fixing the macroturbulence for the stars 
of a given color, without modeling the stars individu- 
ally. If the macroturbulence value was underestimated 
for T > 5800 K, the resulting vsini values (especially 
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Fig. 5.— Rotatio nal velocity histogram for 276 F8-K2 (0.9 < 
M < 1.1 Mq) stars dValenti fc Fi scher 2005). The Sun (I'siniQ = 
1.28 km s"-*) rotates more slowly than 83 ±7% of the stars. There 
is 1 star to the right of the plot with vsini = 36 km s~^. 



of lValenti fc Fischej[2005h . 

The incHnation of the stellar rotational axis to the line 
of sight is usually unknown so the observable is usini. 
Using the solar sp e ctruin reflected by the asteroid Vesta, 
IValenti fc FischeJ (|2005D derived a solar wsinz = 1.63 
km s~^. For the purposes of this analysis we use the 
mean value that would be derived for the Sun, when 
viewed from a random inclination: iisinjQ = 1.63(7r/4) 
km s^i w 1.28 km s"^ 

The Sun rotates more slowly than 83 ± 7% of the stars 
in our Valenti fc Fischer (2005 ) sarn ple (Fig. O. This is 
in agreement with ISoderblomI ( 19831 .11985) who reported 
that the Sun is within 1 standard deviation of stars of its 
mass and age. 

2.6. Galactic Orbital Parameters 

The Galactic velocity components of a star {U,V,W) 
with respect to the local standard of rest (LSR) may 
be used to compute a star's orbit in the Galaxy. How 
typical or atypical is the solar orbit compared to the or- 
bits of other nearby stars in the Galaxy? The orbit may 
be related to habitability because more eccentric orbits 
bring a star closer to the Galactic center where there is 
a larger danger to life from supernovae explosions, cos- 
mic gamma and X-ray radiation a nd any factors associ- 
ated with higher stella r densities (jGonzalez et al.l [20011 : 
iLineweaver et aLll2004D . 

For a standard mo del of the Galactic potential, 
iNordstrom et al.l (120041 ) computed orbital paramters for 
the Sun, and for a large sample ('^ 16700) of A5-K2 stars. 
Their adopted components of the solar velocity relative 
to the local standard of rest were ( U, V, W) = (10.0 ± 
0.4,5 .25 ± 0.62, 7.17 ± 0.38) km s'^ (jPehnen fc BinnevI 
[T998h . 

For each of the 19 87 stars within 40 pc in the 
INordstrom et al.l ()2004f ) catalog, an inner and outer radii 
i?i„in and i?max wcrc computed. This yielded the or- 




Fig. 6. — Top : eccentricity distributio n for the 1,987 stars at 
d < 40 pc from INordstrom et al.l I I2004I '). The Sun has a more 
circular orbit than 93 ± 1% of the A5-K2 stars within 40 pc. Af- 
ter mass, eccentricity is the second most anomalous parameter. 
Bottom panel: Correlation between iilsr a^nd eccentricity for the 
same stars presented in the Bottom: since these properties are 
highly correlated we select only 1 for the analysis. The large grey 
point with error bars represents the median and the 68% widths of 
the two one-dimensional distributions. As in Fig. |4] the contours 
correspond to 38%, 68%, 82% and 95%. 



nents of the solar motion (jPehnen fc BinnevI Il998f l rel- 
ative to the local sta ndard of rest in the Galactic po- 
tential of iFlvnn et"all ([l996). The bottom panel of Fig- 
ure [6] shows the correlation between Galactic orbital ec- 
centricity e and the magnitude of the galactic orbital 
velocities with respect to the local standard of rest: 
ulsr = (U"^ + V'^ + W'^Yl'^ . Eccentricity e and wlsr are 
strongly correlated. We include e, not i^lsKj in the anal- 
ysis since e is less correlated with the maximum height 
above the Galactic plane 2^max, than is wlsr- This is 
shown in Fig.[Tn|in Appendix A. 

The Sun's eccentricity was determined with the same 
relation as the stellar eccentricities. The uncertainty in 
our estimate of solar eccentricity came from propagating 
the uncertainty in the adopted solar motion. We find 
60 = 0.036 ± 0.002 (consistent wi th the Cq = 0.043 ± 
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40 pc (with median eccentricity /Lti/2 = 0.1). This is 
the second most anomalous of the 11 solar properties we 
consider here. 

The frequency of the passage of a star through the thin 
disk could be associated with Galactic gravitational tidal 
perturbations of Oort cloud objects that might increase 
the impact rate on potentially habitable planets. This is 
correlated with the maximum height, Z^ax, to which the 
stars rise above the Galactic plane. Figure [7] shows the 
stellar distribution of Z^ax for the stars shown in Figure 
[SI We find that 59 ± 3% of the A5-K2 stars within 40 
pc of the Sun reach higher above the Galactic plane than 
the Sun does (^max,© = 0.104 ± 0.006 kpc). The solar 
^max,0 was derived by integrating the solar orbit in the 
Galactic potential. The uncertainty on W, produces the 
uncertainty on Zmax,© and hence the ±3% uncertainty on 
59%. Our results for ecce ntricity and ^max a re consistent 
with those obtained using iHogg et all (iOOSi) LSR values: 
(U, V, W) = (10.1 ± 0.5, 4.0 ± 0.8, 6.7 ± 0.2). Using the 
Hogg et al. LSR values, 92 ± 1% of A5-K2 stars within 
40 pc have higher eccentricities than the Sun and 62 ±4% 
of A5-K2 stars within 40 pc have larger Z^ax values. 

How does the Sun's distance from the center of the 
Milky Way compare to the distances of other stars from 
the center of the Milky Way? In Fig. [5] we show the 
distribution of the mean radial distances of stars from 
the Galactic center, bas ed on the star count model of 
iBahcall fc Soneiral (|198Cl( ). To represent the entire Galac- 
tic stellar population we include the disk (thin + thick) 
and spheroidal (bulge + halo) components. Using the 
curre nt Solar distance f rom t he center (i?o = 7.62 ± 0.32 
kpc, lEisenhauer et al.l I2005D and a disk scale length 
h = 3.0 ± 0.4 kpc (iGould et aI]|1996D . we estimate that 



the Sun lies farther from the Galactic center than 72l5% 
of the stars in the Galaxy. The uncertainty on the re- 
sult comes from the 68% bounds of the total distribution, 
which come from the scale length uncertainty (±0.4 kpc). 




Fig. 7. — The distri bution of maximu m heights above the Galac- 
tic plane for the ^ordstrom et al.l II2004I) sample. 59 ±3% of nearby 
A5— K2 stars (dmax =40 pc) reach higher above the Galactic plane 
than the Sun reaches. There are 22 stars evenly distributed over 
^max between 1.5 and 9.6 kpc. Their exclusion from the compari- 
son reduces the 59% result by less than 1%. 




[kpc] 



Fig. 8. — Mean stellar galactocentric radius distribution 
dNi, / dRQ^\. The solid curve represents the sum of the disk (dashed 
line) and spheroidal (dotted line) stellar components. The 68% 
uncertainty of the total distribution is shown by the cross-hatched 
area, 

the stars in the Galaxy. 



The Sun is farther from the Galactic center than 72^g% of 



2.7. Host Galaxy Mass 

The mass of a star's host galaxy may be correlated 
with parameters that have an influence on habitability. 
For example, galaxy mass affects the overall metallicity 
distribution that a star would find around itself — an 
effect that would not show up in Fig. [31 which only shows 
the local metallicity distribution. 

The Milky Way is more massive than ~ 99% of all 
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log,o Host Galaxy Stellar Mass [Mq] 

Fig. 9. — Fraction of all stars that live in galaxies of a given 
mass, dNi,/dM (solid curve). The mass of the Sun's galaxy is 
indicated by the "0". This distribution represents the amount 
of stellar mass contributed by galaxies of a given mass. Approx- 
imately 77^^^^% of stars live in galaxies less massive than ours. 
The cross-hatched band shows the la uncertainty associated with 
the uncertainty in the two Sc hechte r function parameters, a and 
L* llLwcdav 2000; Schcchte^ [T976h . The dashed line shows the 
unweighted luminosity function (the number of galaxies per lumi- 
nosity interval dNg^ia/dM) according to which the Milky Way is 
more massive than ~ 99% of galaxies. 



galaxy, and the beha viour of th e low-mass end of the 
galaxy mass function (|Silkll2007D . We are referring here 
to the stellar mass, not the total baryonic mass or the to- 
tal mass. Despite the Milky Way's large mass compared 
to other galaxies, if most stars in the Universe resided 
in even more massive galaxies, the Milky Way would be 
a rather low mass galaxy for a star to belong to. To 
estimate the fraction of all stars in galaxies of a given 
mass, we first estimate the distribution of galax y masses 
by taking the K-hand luminosity function of iLovedavl 
(pOOO.) (iiT-band most closely reflects stellar mass since it 
is less sensitive than other bands to differences in stellar 
populations) and weighting it by luminosity. We convert 
this to stellar mas s assuming a co nstant stellar-mass-to- 
Ught ratio of 0.5 ( Bell fc de Jo^ [20011. This function, 
plotted in Fig. ^ shows the amount of stellar mass con- 
tributed by galaxies of a given mass — or assuming iden- 
tical stellar populations — the fraction of stars residing 
in galaxies of a given stellar mass. 

We estimate the _?C-band luminosity of the Milky 
Way by converting the published y-band magnitude of 
iCourteau &: van den Berghl (|1999l ) to the iiT-band assum- 
ing the mean color of an Sbc spira l galaxy from the 
2 MASS Large Galaxy Atlas (iJarrett e t al. 2003) and ap- 
plying the color conversion from (jPriver et al.lll994ll . We 
then convert this to stellar mass using the same stellar- 
mass-to-light ratio used above, i.e., 0.5. In this way we 
estimate the stellar-mass content of the M ilky Way to 
be io io-55±o.i6 ^ s.etil X IQI^Mq (see also lFlvnn etall 
l2006f) . Comparing this to the stellar masses of other 
ffa.la.xies (Fiff. 151). we find that 77~^}]% of stars reside in 



2.8. Host Group Mass 

The mass of a star's host galactic group or galactic 
cluster may be correlated with parameters that have 
an influence on habitability. For example, group mass 
is correlated with the density of the galactic environ- 
ment (number of galaxies per Mpc'^) that could, like 
galactocentric radius, be associated with the dangers of 
high stellar densities: "The presence of a giant ellipti- 
cal at a distance of 50 kpc would have disrupted the 
Milky Way Galaxy, so that human beings (and hence 
astronome rs) probably wovild not have come into ex- 
istence." (jvan den Bergh 200 0). Our Local Grou p of 
galaxies seems rather typical (jvan den Berghl [200(1 ) but 
we would like to quantify this. Proceeding similarly 
to our analysis of galaxy mass in Sect. 12.71 we ask, 
what fraction of stars live in galactic groups less massive 
than our Local Group? Figure [TO] shows the luminosity- 
weighted (i.e. stellar mass-weighted) number density 
of galactic groups. The number distribution and lumi- 
nosity distribution of galactic groups is taken from the 
Two-degree Field Galaxy Redshift Survey Percolation- 
Inferred Galaxy Group (2PIGG) catalogue (Eke et aT] 
I2004D . It spans the range from weak groups to rich galaxy 
clusters. 

We estimated the stellar masses of the 2PIGG groups 
and Local Group galaxies (jCourteau fc van den Bergh[ 
|1999() by converting from the B-band a ssuming a con- 
stant stellar-mass-to-light ratio of 1.5 (jBell fc de Jongj 
[200l . This gives an estimated stellar mass of the lo- 
cal group of I0i0-9i±" 07 ^ 8.lt\i X 10^°Mq. Figure 
[TOl indicates that our Local Group is a typical galactic 
grouping for a star to be part of. Approximately 58 ±5% 
of stars live in galactic groups more massive than our 
Local Group. With respect to the mass of its galaxy and 
the mass of its galactic group, the Sun is a fairly typical 
star in the Universe. 

3. JOINT ANALYSIS OF 11 SOLAR PROPERTIES 

3.1. Solar x^- Analysis 

We would like to know whether the solar properties, 
taken as a group, are consistent with noise, i.e., are they 
consistent with the values of a star selected at random 
from our stellar distributions. We take a approach to 
answering this question. First we estimate the solar Xg, 
by adding in quadrature, for all 11 properties, the dif- 
ferences between the solar values and the median stellar 
values. We find: 

N=ll I x2 
2 _ l^0.« ~ ^1/2,^ _ ^ OO+0.08 

Xq 



i=l 



7.: 



-0.30 



(1) 



where i is the property index, N = 11 is the number of 
properties we are considering, /ii/2,i is the median of the 
i*'' stellar distribution and (Je8,i is the difference between 
the median and the upper or lower 68% zone, depending 
on whether the solar value XQ^i is above or below the 
median. The uncertainty on Xq is obtained using the 
uncertainties of XQ^i. 

Equation ([T]) can be improved on by taking into ac- 
count i) the non-Gaussian shapes of the stellar distribu- 
tions and a) the larger uncertainties of the medians of 
smaller samples (our smallest sampl e is ^100 s tars). 
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Fig. 10. — The dashed histogram shows the luminosity function 
of galactic groups (number of groups per interval of B-band lumi- 
nosity). The solid histogram shows the luminosity- weighted group 
luminosity function (approximately the fraction of stars which in- 
habit a group of given stellar mass). The horizontal axis has been 
converted to stellar rnass assuming a constant B-band stellar-mass- 
to-light ratio of 1.5 (Bell & de Jong 2 00l|). The "0" shows t he esti- 
mated mass of the Local Group tCourteau &: van den Berghll 19991) 
and lies just below the median (vertical grey line). 



more accurate estimate of Xq- Because the bootstrap is 
a non-parametric method, the distributions need not be 
Gaussian. 

We obtain x| = 8.39 ± 0.96. Figure [U] shows the 
resulting solar distribution. The median of this dis- 
tribution is our adopted solar value. Dividing our 
adopted solar by the number of degrees of freedom 
gives our adopted reduced solar value: 



x|/ll = 0.76 ±0.09 



(2) 



The standard conversion of this into a probability of 
finding a lower value (assuming normally distributed 
independent variables) yields: 



P(< x| = 8.39|A^ = 11) = 0.32 ± 0.09. 



(3) 



3.2. Estimate of P{< x|) 

To quantify how typical the Sun is with respect to our 
11 properties, we compare the solar Xq{— 8.39) to the 
distribution of x^ values obtained from the other stars in 
the samples. 

We perform a Monte Carlo simulation 

([Metropolis fc UlamI Il949( ) to calculate an estimate 
of each star's x^ value (x*)- The histogram shown 
in Figure [12] is the resulting Monte Carlo stellar x^ 
distribution. Three standard x^ distributions have 
been over-plotted for comparison {N = 10, 11, 12). The 
probability of finding a star with x^ lower than or equal 
to solar is: 

aic(<X0 =8.39|Ar= 11) = 0.29 ±0.11 (4) 
The Monte Carlo x^ distribution has a similar shape 




Solar X 

Fig. 11. — Bootstrapped solar distribution, 
the distribution (white "0") is Xq = 8-39 ± 0.96. 

compared to the solar Xq value from Eq. [T] 7.S 
over-plotted (grey "0" on dotted line). 



The median of 
This should be 



0.29 ~ P{< x^) = 0.32 (Eqs.[3]and|4l). The more appro- 
priate Monte Carlo distribution has a longer tail, pro- 
duced by the longer super-Gaussian tails of the stellar 
distributions. 

Table [2] summarizes our analysis for the solar Xq values 
and the probabilities P{< x©)- Our simple x© = 7.88 
estimate increased to 8.39 and the uncertainty increased 
by a factor of ~ 3 after non-Gaussian and sample size 
effects were included as additional sources of uncertainty. 
Our improved analysis yields Pmc(< Xq): with a longer 
tail and brings the probability down from 0.32 ± 0.09 
to 0.29 ± 0.11. If this value were close to 1, almost all 
other stars would have lower x^ values and we would have 
good reason to suspect that the Sun is not a typical star. 
However, this preliminary low value of 0.29 indicates that 
if a star is chosen at random, the probability that it will 
be more typical (~ have a lower x^ value) than the Sun 
(with respect to the 11 properties analysed here), is only 
29 ± 11%. The details of our improved estimates of x© 
and P(< X©) can be found in the Appendix [6l 

4. RESULTS 

Figure [T3] shows four different representations of our 
results. Figure [T3l A compares the solar values to each 
stellar distribution's median and 68% and 95% zones. 
The Sun lies beyond the 68% zone for three properties: 
mass (95%), eccentricity (93%) and rotational velocity 
(88%). No solar property lies beyond the 95% zone. The 
histogram in FigurjTSB is the distribution of solar values 
in units of standard deviations: 



Ml/2,j 



(5) 



For each stellar property i, the Sun has a larger value 
than ni% of the stars. If the Sun were a randomly se- 
lected star, we would expect the percentages ni% to be 
scattered rouffhlv evenlv between 0% and 100%. 
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TABLE 2 

Summary of and P(< x%) results. 



Analysis 


X% 




P{<X%\N = 11) 


^'mc(<x|,I^' = ii) 


simple 
improved 


7.88lO:»g (Eq.[TJ 
8.39 ±0.96 


n 72+0 01 

'^-0.03 
0.76 ±0.09 (Eq.EJl 


0.28l0 »i (Eq.0 
0.32 ±0.09 (Eq.|3]l 


0.29 ±0.11 (Eq.mi 
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Fig. 12. — Stellar x'^ distribution from our Monte Carlo simula- 
tion. Pmc(< Xq = 8-39) = 0.29 ± 0.11 (represented by the grey 
shade) is calculated integrating from = to = Xq- For 
comparison, three x^ distribution-curves are over-plotted with 10, 
11, and 12 degrees of freedom. The standard probability from the 
iV = 11 curve yields: P(< x| = 8.39|Ar = 11) = 0.32 ± 0.09. 
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Fig. 13. — Various representations of our main results. A: Solar values of 11 properties compared to the distribution for each property 
Each distribution's median value is indicated by a small filled circle. The dark and light grey shades represent the 68% and 95% zones 
respectively. B: Histogram of the number of properties as a function of the number of standard deviations the solar value is from the 
median of that property. The grey curve is a Gaussian probability distribution normalised to 11 parameters. C: Percentage 71^% of stars 
with sub-solar values as a function of property. The average signal expected from a random star is shown by the solid line (see Sec.|4]l. D: 
Percentage % of stars with sub-solar values as a function of the number of standard deviations the solar value is from the median of that 
property. The solid curve is a cumulative Gaussian distribution — if every sample were a Gaussian distribution, every solar dot would sit 
exactly on the line. Just as in (C), the dashed lines encompass the 68% and 95% zones. Similar to the results from Figure [T2l these four 
panels indicate that the Sun is a typical star. 
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When the ni% values are lined up in decreasing order 
fFig.fOlC). we expect them to be near the line given by: 



^i, expected^ 



(z-1/2) 



N 



X 100% 



(6) 



and plotted in Figure fTSl C . Any anomalies would show up 
as points '0' significantly distant from the line. 

FigurdTSD) compares the percentages ni% of stars 
having sub-solar values (shown in Figurc[l3P) with the 
solar values expressed in units of standard deviations 
from each distribution's median (shown in Fig. I13B). 
If the stellar distributions were perfect Gaussians, the 
translation from Zi to rii would be given by the cumula- 
tive Gaussian distribution (Fig. I13P. solid curve). That 
the points lie along this line demonstrates that the ap- 
proximation of our distributions as Gaussians is reason- 
able. 

Table [3] lists percentages ni% of stars for each property 
(as shown in Fig. [T5|) . In the lower half of the table we 
list properties not included in this analysis because of 
correlations with properties that are included. 

Individual stellar uncertainties make the observed 
characteristic widths (cres, Table [T]) larger than the 
widths of the intrinsic distributions. This broadening 
effect makes the Sun appear more typical than it really 
is when ergs and the individual stellar uncertainties (cr^) 
are of similar size and the individual stellar uncertainties 
are much larger than the solar uncertainty (ctq). We es- 
timate that our results are not significantly affected by 
this broadening effect. 

Our resulting probability of finding a star with a 
lower or equal to the solar value of 29 ± 11% (Eq. [4|), is 
consistent with the probability we would obtain if stellar 
multiplicity were included in our study. Using the vol- 
ume limited sample used for stellar mass in Section [2T] 
(125 A1-M7 stars within 7.1 pc) the probability that 
a randomly selected star will be single is 52.8 ± 4.5%, 
which means that ~ half of stars are single while the 
other half have one or more companions. Including this 
in our bootstrap analysis and Monte Carlo simulations 
(see Appendix [6]) marginally increases the probability in 
Eq. H to 33 ±11%. If the multiplicity data for 246 G 
dwarfs from 'D uquennov fc Mavorj (|1991f ) is used instead 
— the probability that a randomly selected G dwarf will 
be single is 37.8 ± 2.9% — then the probabihty in Eq. H 
would increase to 34 ± 11%. The inclusion of stellar mul- 
tiplicity marginally incr eases our repo r ted pr obability. 

In Figures 6 and 7 of iRadick et al.l (|1998D . the Sun's 
short-term variability as a function of average chromo- 
spheric activity, appears ~ la low, co mpared to a dis- 
tribu tio n of 35 F3-K7 S un-lik e stars (|Lockwood et al.l 
Il997( ). iLockwood et al.l ()2007l ) suggest that the Sun's 
small total irradiance variation compared to stars with 
similar mean chromospheric activity, may be due to their 
limited sample and the lack of solar observations out of 
the Sun's equatorial plane. We do not include short or 
long term variability (chromospheric or photometric) in 
Table [3] because of the small size of the ILockwood et al.l 
([200/1) sample. We also do not include the chromospheric 
index i?HK (^^^ Tabled bottom panel) as one of our 11 
properties because of its correlation with the chromo- 
spheric ages of our sample. 



5. DISCUSSION AND INTERPRETATION 

The probability Pmc(< x|) = 0.29 ± 0.11 classifies 
the Sun as a typical star. How robust is this result? 
The probability of finding a star with a lower than or 
equal to Xq: depends on the properties selected for the 
analysis (see problem Hi of Section I). For example, if we 
had chosen to consider only mass and eccentricity data, 
this analysis would yield Pucix^ < Xq) = 0-94 ± 0.4, 
i.e., the Sun would appear mildly (~ 2a) anomalous. 
If on the other hand, we had chosen to remove mass 
and eccentricity from the analysis, we would obtain 
Pucix^ ^ Xq) = 0-07 i 0.04, which is anomalously low. 
The most common cause of such a result is the over- 
estimation of error bars. The next most common cause is 
t he preselecti o n of p roperties known to have Tii% ^ 50%. 

iGustafssonl ((1998) discussed the atypically large so- 
lar mass, and proposed an anthropic explanation — the 
Sun's high mass is probably related to our own existence. 
He suggested that the solar mass could hardly have been 
greater than ~ 1.3 since the main sequence l ifetime of 
a 1.3 Mq star is 5 billion years l|Clavtonl ll983). He also 
discussed how the dependence of the width of the circum- 
stellar habitable zone on the host star's mass probably 
favours host stars within the mass range 0.8-1.3 M0. 

Our property selection criteria is to have the largest 
number of maximally independent properties that have 
a plausible correlation with habitability and, ones for 
which a representative stellar sample could be assem- 
bled. Our joint analysis does not weight any parameter 
more heavily than any other. If the only properties as- 
sociated with habitability are mass and eccentricity then 
we have diluted a ~ 2(t signal that would be consistent 
with Gustafsson's proposed anthropic explanation. 

Our analysis points in another direction. If mass and 
eccentricity were the only properties associated with hab- 
itability, then the solar values for the remaining nine 
properties would be consistent with noise. However, a 
joint analysis of just the remaining nine properties pro- 
duces a x© 9 = 3.6 ±0.4 and the anomalously low proba- 
bility: P(<' x|,9) = 0.07 ±0.04, which suggests that the 
nine properties are unlikely to be the properties of a star 
selected at random with respect to these properties. 

The x^ fit of the 11 points in Panel C of Fig. 13 to the 
diagonal line yields a fit that is substantially better then 
the fit of the remaining nine properties to Eq. 7 with 
= 9. In other words, the joint analyis suggests that 
although mass and eccentricity are the most anomalous 
solar properties, it is unlikely that they are associated 
with habitability, because without them, it is unlikely 
that the remaining solar properties are just noise. Thus, 
the Sun, despite its mildly (~ 2a) anomalous mass and 
eccentricity, can be considered a typical, randomly se- 
lected star. 

There may be stellar properties crucial for life that 
were not tested here. If we have left out the most impor- 
tant properties, with respect to which the Sun is atypical, 
then our Sun-is-typical conclusion will not be valid. If we 
have sampled all properties associated with habitability, 
our Sun-is-typical result suggests that there are no spe- 
cial requirements on a star for it to be able to host a 
planet with life. 
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TABLE 3 

Summary of How the Sun Compares to Other Stars (see Fig. [TSl l 



Parameter 


Figure 




Level of Anomaly 


Mass 


m 




95 ± 2% 


of nearby stars are less massive than the Sun. 


Age 


E 




53 ± 2% 


of stars in the thin disk of the Galaxy are older than the Sun. 


[Fe/H] 


11 




65 ± 2% 


of nearby stars are more iron-poor than the Sun. 


[C/O] 


m 


81 ± 3% 


of nearby stars have a higher C/O ratio than the Sun. 


[Mg/Si] 


HP 


66 ± 3% 


of nearby stars have a higher Mg/Si ratio than the Sun. 


V sin i 


E 




83 ± 7% 


of nearby Sun-like-mass stars rotate faster than the Sun. 


e 


E 




93 ± 1% 


of nearby stars have larger galactic orbital eccentricities than the Sun. 


^max 


E 




59 ± 3% 


of nearby stars reach farther from the Galactic plane than the Sun. 


R-Gal 


E 




72j;^% 


of stars in the Galaxy are closer to the galactic center than the Sun. 


Mgai 


i 




77ll^% 


of stars in the Universe are in galaxies less massive than the Milky Way. 


group 


[El 


58 ± 5% 


of stars in the Universe are in groups more massive than the local group. 



Properties Not Included in the Analysis Because They Are Correlated with the Selected 11 Parameters 



Mass: IMFstoiiar 


m 


94 ± 2% 


of nearby stars are less massive than the Sun. 


Age: SFRcosmic 


m 


86 ± 5% 


of stars in the Universe are older than the Sun. 


Age- 




55 ± 2% 


of nearby Sun-like-mass stars are older than the Sun. 


[Fe/K]^ 




56 ± 5% 


of nearby stars are more iron-poor than the Sun. 


V sin i'^ 




92 ± 5% 


of nearby Sun-like-mass stars rotate faster than the Sun. 


logR' d 

[0/Fe1 




51 ± 2% 


of nearby FGKM stars are more chromospherically active. 




75 ± 3% 


of nearby stars have a lower O/Fe ratio than the Sun. 


-^min 




91 ± 1% 


of nearby stars get closer to the Galactic center. 


fLSR 




93 ± 1% 


of nearby stars have smaller velocity with respect to the LSR. 


\u\ 




75 ± 1% 


of nearby stars have larger absolute radial velocity. 


\v\ 




82 ± 1% 


of nearby stars have larger absolute tangential velocity. 


\w\ 




58 ± 1% 


of nearby stars have larger absolute vortical velocity. 



- 1126 stars (A5-K2) from 'Nordstrom et al.| ll200i '). 
^ 91 stars (GK) from Fivata ct al. (1997). 
" 590 stars (F8-K2) from Nordstrom ct al? f2Q0D. 
866 stars (FGKM) from Wrigh t _c,t a l, (2004.) . 

6. CONCLUSIONS 

We have compared the Sun to representative steUar 
samples for 11 properties. Our main results are as fol- 
lows: 

• Stellar mass and Galactic orbital eccentricity are 
the most anomalous properties. The Sun is more 
massive than 95 ± 2% of nearby stars and has a 
Galactic orbital eccentricity lower than 93 ± 1% 
FGK stars within 40 pc. 

• Our joint bootstrap analysis yields a solar Xq — 
8.39 ± 0.96 and a solar reduced X^/ H = 0.76 ± 
0.09. The probability of finding a star with a 
lower than or equal to solar i-Mc(< Xq — 8.39 ± 
0.96) = 0.29 ±0.11. 

To our knowledge, this is the most comprehensive and 
quantitative comparison of the Sun with other stars. We 
find that taking all 11 properties together, the Sun is a 
typical star . This finding is largely in agreement with 
iGustafssonI ()1998l ). however our results undermine the 
proposition that an anthropic explanation is needed for 
the comparatively large mass of the Sun. 

Further work could encompass the inclusion of other 
properties potentially associated with habitability. An- 
other improvement would come when larger stellar sam- 
ples become available for which all properties could be 
derived, instead of using different samples for different 
properties as was done here. In addition, research in the 
molecular evolution that led to the origin of life may, in 



on Earth, orbiting the Sun. 
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APPENDIX A 
PROPERTY CORRELATIONS 

The formalism and the use of the ^^-distribution 
to obtain P{< Xq\N), — improved using Monte-Carlo 
simulations in Section [3?2l to obtain Pmc(< Xq) — as- 
sumes that each parameter is independent of the others. 
In selecting our 11 properties we have selected proper- 
ties which are maximally independent based on plotting 
property 1 vs property 2 for the same stars. Wc show 
seven such plots in this Appendix. 

If there are correlations between the analysed proper- 
ties, then the number of degrees of freedom N could drop 
from 11 to 10.5 (see Fig. 12). Some properties have 
been excluded from the analysis due to a correlation with 
another property in the analysis. 

Al. ELEMENTAL RATIOS 
In Figure [14] we show the distribution for carbon to 
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Fig. 14. — Carbon to oxygen ratio [C/O] versus magnesium to 
silicon ratio [Mg/Si] of 176 FG stars with abundances for these ele- 
ments ( ReddY et al.,:2003). In Figs.|4p and|4p we showed that the 
[C/O] and [Mg/Si] distributions are largely independent of [Fe/H]. 
Here we show that these distributions are also largely independent 
of ea ch other. Note that in this comparison we only use the data 
from IReddv et al.l 120031) , since it is the largest available sample 
with C, O, Mg and Si abundances. 

A2. MASS, AGE, AND ROTATIONAL VELOCITY 

In Figure [15] we show four correlation plots for mass, 
chromospheric age, rotational v elocity and f sin t. We 
use the stars common t o both Wrig ht et al.l ()2004f ) and 
iValenti fc Fischeil (|2005l) for which these observables are 
available. 

A3. GALACTIC ORBITAL PARAMETERS 

The Galactic orbital eccentricity (e) and the magnitude 
of the galactic orbital velocities with respect to the local 
standard of rest (i^lsr) are strongly correlated (see Fig. [6] 
in Sec. 12. 6p . We selected e instead of ulsr because of 
its near independence of the maximum height above the 
galactic plane (.^max); see Figure fTBI 

APPENDIX B 
IMPROVED ESTIMATES OF Xq AND P(< Xq) 

In Section 13.21 with 11 degrees of freedom, the re- 
duced from EquationHis / 11 = Since 
X0 / 11 < 1, the Sun's properties are consistent with the 
Sun being a randomly selected star. 

To improve on this preliminary analysis (but with a 
similar conclusion), as me ntioned in Se ction we em- 
ploy a bootstrap analysis (jEfronl fl979l ) to randomly re- 
sample data (with replacement) and derive a more ac- 
curate estimate of x©- Because the bootstrap is a non- 
parametric method, the distributions need not be Gaus- 
sian. 

For every iteration, each parameter's stellar distribu- 
tion is randomly resampled and a x© value is calculated 
using Eq. ([T]). The uncertainties Gq^i of the solar values 
X0_i are also included in the bootstrap method: for every 
iteration, the Solar value for each parameter is replaced 
in Eq. ([T]) by a randomly selected value from a normal 



although the resulting distribution varies very little once 
the number of iterations reaches ^ 10,000. 

The median of this distribution and the error on the 
median yields our improved value for the reduced x© 
(Fig. [TT|) . The uncertainty of the median of each re- 
sampled distribution varies inversely proportionally to 
the square root of the number of stars in the distribu- 
tion, A/i]^/2,i oc 1/ yjN\~i. In other words, median values 
are less certain for smaller samples and this uncertainty 
is included in our improved estimate of x© ■, and its un- 
certainty. 

We find the probability of finding a star with a x* value 
lower than the solar x©, for N = 11 degrees of freedom 
in the standard way (jPress et al.lll99^ and obtain: 



BlH|ll) = 0.28t° 



-fO.Ol 
03 



(7) 



To improve our estimate of the probability of finding 
a star with lower x^ valu e than the Sun, we per form 
a Monte Carlo simulation (jMetropolis fc Ulamlll94 9) to 
calculate an estimate of each star's value (x*)- For 
every iteration, we randomly select a star from each stel- 
lar distribution. We then calculate its x* value by re- 
placing the solar value x©^i with that star's value Xi,^i 
in Eq. ([1]). This process was repeated 100,000 times to 
create our Monte Carloed stellar x^ distribution. Stars 
were randomly selected with replacement, thus the sim- 
ulated x^ distribution accounts for small number statis- 
tics and non-Gaussian distributions. The probability of 
finding a star with x^ lower than or equal to solar is 
Pmc = 0.29 ±0.11. 

The results of our analysis for the Solar x© values and 
the probabilities P(< x©) are summarized in Table [2] 

Bl. ADDITION OF A DISCRETE PARAMETER 

In Section[5]we discuss the addition of stellar multiplic- 
ity to our analysis. Since stellar multiplicity cannot easily 
be approximated by a one-sided Gaussian (particularly 
because the Sun is on the edge of the distribution, i.e., 
it is of multiplicity one), we modified our Monte Carlo 
procedure to include this discrete parameter. The likeli- 
hood of observing a particular x^ for the 11 parameters 
is 

expf-iVxA. (8) 



We take the probability p(l) of a star being a single 
star, to be 53.8 ± 4.5%, obtained from our sample of 
nearby stars (Sec. 12. l| ). The likelihood L of observing a 
particular x^ and p{\) is the product 



(9) 



Taking logarithms we can then compute the distribution 
of the statistic 5, where 



(10) 



/-1 1 ci-J--**! V\1 t-i-ti- 
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Fig. 15. — Correlation plots between various properties. For all four panels we use the stars common to both I Wright et al.l (120041 ) and 
IValenti fc Fischeri l|2005h . P anel (A ) I: ma ss vs rotational velocity vsini for 713 FGK stars. This panel shows the degree of correlation 
See IGravl II2005I ) for a stronger correlation between these two variables when a larger mass range and more active 
To minimize the effect of this correlation on our analysis, we restrict the range of mass in Fig. 5 to 0.9 to 

The lack of correlation between chromospheric determined ages 
and rotational velocities is shown. Panel (C): no strong correlation between mass and chromospheric age for 639 FGK stars. Panel (D): 
the ages of 637 stars determined by the chromospheric method versus their ages from the isochrone method. 



between mass and i;sini. 
stars are kept in the sample. 
I.IMq. Panel (B): chromospheric age versus usini for 641 FGK stars. 
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